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FOREWORD 

The  following  report  will  develop  a  complete  set  of  field  equations  for  the  dynamic  motions 
of  a  towed  cable  system.  These  equations,  when  solved,  constitute  a  dynamic  simulation  of  a 
towed  cable  system. 

During  gentle  maneuvers,  the  important  forces  in  the  governing  force  balances  on  a  cable 
system  are  tension  and  hydrodynamic  drag.  In  more  severe  maneuvers  involving  highly  dynamic 
situations,  however,  the  effects  of  shear,  bending,  torsion,  extension,  rotary  inertia,  and  inertia  can 
be  inqxntant,  and  sometimes  even  mathematically  necessary. 

Many  dynamic  models  for  simulating  the  modons  of  a  towed  cable  system  already  exist 
Many  of  these  models,  however,  by  ignoring  the  above-mentioned  effects,  have  lestrictive 
assumptions  built  into  them.  The  model  represented  by  the  equations  develqped  here  does  not 
have  these  restrictions.  Thus,  the  equations  can  be  used  either  for  more  accurate  simulations  of 
highly  dynamic  situations  or  simply  for  studying  the  range  of  validity  of  simpler  models. 


vii/viii 
Reverse  Blank 
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A  COMPLETE  SET  OF  FIELD  EQl  TIONS 
FOR  THE  DYNAMIC  SIMULATION 
OF  A  TOWED  CABLE  SYSTEM 

THE  COORDINATE  SYSTEMS 

The  total  field  equations  describing  the  dynamic  behavior  of  a  towed  system  involve  four 
coordinate  systems;  1)  a  fixed  inertial  coordinate  system  in  which  to  write  the  translational 
dynamic  equations,  2)  a  body-fixed/principai-axis  coordinate  system  in  which  to  write  the 
rotational  dynamic  equations.  3)  a  flow-based  coordinate  system  in  which  to  write  the 
hydrodynamic  loading,  and  4)  a  coordinate  system  based  on  the  local  cable  curvature  in  whicli  to 
express  the  bending  and  twisting  moments. 


THE  BODY-FIXED  COORDINATE  SYSTEM  (i,2,3) 

Consider  two  orthonormal  coordinate  systems;  a  fixed  inertial  one  (x,  y,  z) ,  and  a  body- 
fixed  one  (i,  2,  §].  (See  figure  1.)  Three  successive  rotations,  through  the  Euler  angles  (p,  0,  and 
Y,  can  be  performed  to  transform  (x,  y,  z)  to  ^1, 2,  sj. 
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A 

z 


Figure  1.  The  Inertial  and  Body-Fixed  Systems 

The  first  rotation,  through  the  angle  <p  about  the  z  axis,  iransforms  the  (x,  y,  z)  system  to 
an  (x',  y',  z')  system  where,  of  course,  z  coincides  with  z'.  The  second  rotation,  through  the 
angle  0  about  the  x'  axis,  transforms  the  (x',  y',  z')  system  to  an  (x",  y",  z")  system  where,  of 
course,  x''  coincides  with  x'.  Rnally,  the  third  rotation,  through  the  angle  y  about  the  z'  axis, 
transforms  the  (x'',  y*  z")  system  to  our  body-fixed  system  (i,  2, 3).  Figures  2  through  4  show 
these  rotations,  as  well  as  the  resulting  rotational  matrices  Tl^,  Hq,  and  Tty. 


2 


*<> 


TR  10.231 


cosq> 

-sinq) 

0 


sin<p  0 

COS9  0 

0  I 


“'(x'y'z') 


=  y[xy».) 


Figure  2.  The  ip  Rotation 


0  0  ■ 
cos8  sin  6 
-sin0  COS0 


Figure  3.  The  $  Rotation 


cosy 

“Siny 

0 


siny 

cosy 

0 


0 

0 

1 


V(123)  =  Tty 


Figure  4.  The  Y  Rotation 


3 


TR  10^31 
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Thus,  the  inertial  and  body-Fued  cocnidinate  systems  are  related  as  follows, 


^(123)  =  ^  %y^) 

(1) 

rt  =  fly  n#  fly 

(2) 

‘  cns 9 cos y -sin 9 cos 6 sin 

sin  9  cos  y -f  cos  9  cos  6  sin  y 

sin  0  sin  y' 

fl= 

-C0S9  sLny  -  sin9  cosd  cosy 

-sin  9  sin  y  -f  cos  9  cos8  cos  y 

sinOcosy  (3) 

sinysinO 

-cosysinO 

COS0 

xi  y*  il 

n»  x-2  y  2  z-2  (4) 

xS  y-3  z-3 

■  • 

and,  in  equations  (3)  and  (4),  we  have  spelled  out  the  components  of  the  !,  2,  and  3  unit  vectors 
in  the  (x,  y,  z)  system.  However,  because  we  will  also  be  needing  the  components  of  !\  2',  and 
3'  in  the  (x,  y,  z)  system,  now  is  a  good  time  to  write  them  out,  too. 

COMPUTATION  OF  i',  2%  AND  3' 

The  derivatives  of  the  cable-fixed  unit  vectors  with  respect  to  unstretched  cable  arc  length 
will  be  derived  here  (they  will  be  needed  in  the  next  section).  By  the  chain  rule  we  can  write 


il 

dl 

B  — 

di 

de 

♦  il 

dy 

dio 

dy 

Bio 

de 

Bio 

dy 

Bio 

B2  ^  92 

iSL  + 

B2 

de 

-  ii 

dy 

BiQ 

dy 

3^0 

de 

Bio 

dy 

Bio 

il  a  ii 

is.  + 

d3 

de 

dy 

3^0 

dy 

3^0 

de 

dfo 

ay 

Bio 

4 
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Thus,  finding  the  derivatives  of  the  (i,  2,  §)  unit  vectors  with  respect  to  Iq  ewnes  down  to  finding 
their  derivatives  with  respect  to  the  three  Euler  angles,  <p,  6,  and  y. 


I 


From  figure  4  we  can  write  from  inspection 


From  figure  3  we  can  write  from  inspection 


dx^ 

n 

80 

u 

^  - 

80 

A 

82" 

80 

-y  • 

(6) 


I 


(7) 


Using  the  rotation  matrix  from  figure  4,  we  can  rewrite  each  double-primed 
coordinate  system  vector  in  equation  (7)  in  terms  of  the  body-fixed  system  unit  vectors.  This 
gives 


’cosy 

-siny  O' 

'8i/80‘ 

■  0 

siny 

cosy  0 

di/de 

= 

A 

3 

0 

0 

83/80 

-sinyl  -cosy2 

I 

I 

I 


which,  when  solved,  yields 
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de 

de 

de 


sin\jr3 
cosy  3 

A  A 

-sinv  1  -  cosy  2, 


(9) 


From  flgure  2,  we  can  write  inspection 


dtp 

dtp 

dsp 


(10) 


Using  the  Hy  and  matrices  from  figures  3  and  4,  we  can  rewrite  each  single-primed 
unit  vector  in  equation  (id)  in  terms  of  the  body-fixed  system  unit  vectors.  This  gives 


cosy 

-siny  0 

dl/dtp 

cose  sin  y  i -f  cose  cosy  2  -  sin  9  3 

cosOsiny 

cosGcosy  -sine 

di/dy 

= 

-cosy  i  + siny  2 

, sine  siny 

^inecosy  cose. 

ds/dtp 

0 

which,  when  solved,  yields 


A 

dtp 

dtp 

ii 

d9 


cose  2  -  sine  cosy  3 
-cosei-fsinesiny  3 
sine  cosy  1  -  sine  siny  2 . 


(12) 
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Finally,  substituting  equations  (6),  (9),  and  (12)  into  equation  (5)  gives  us  the  derivatives  we 
were  after 


il 

dio 


0 


■|  <p'cos0  +  v' 
[-<p'sin0  cosv  +  0'sinv. 


123 


il 

3^0 


-<p' COS0  -  v' 

0 

<p'  sin  0  sin  y  +  0'  cosyj 


123 


il 

9^0 


9' sin0  cosy  -  0' sin y  1 
•  -<p'  sin0  siny  -  0'  cosy  > 


(13) 


THE  FLOW.BASED  COORDINATE  SYSTEM  (n,b,t) 

The  hydrodynamic  forces  on  the  cable  will  be  written  most  easily  in  a  coordinate  system 
using  the  local  flow,  U ,  and  the  cable  to  define  the  system.  The  local  flow  U  is  the  cable  velocity 
through  the  water  and  is  thus  the  caLIe  velocity  u  less  any  currents;  U  =  u  -  U(. .  Figure  5 
defines  such  a  system. 
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A 

n 


A 

b 


Figure  5,  The  Flow-Based  Coordinate  System 

In  this  coordinate  system,  the  tangent  vector  t  is  defined  as  3,  the  same  tangent  vector  used 
in  the  body-fixed  coordinate  system.  Next,  using  t  and  tJ  to  form  a  cable/flow  plane,  the 
k!ri:jimal  b  is  defined  normal  to  the  t/U  plane.  Finally,  using  t  and  b,  the  normal  vectcx:  n  is 
formed,  Crrnpleting  the  CKthonormal  set  ^h,  b,  tj.  Thus,  performing  the  algebra  indicated  in 
figure  3,  and  using  the  results  of  equations  (3)  and  (4),  we  can  write; 
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1 

|U|sina 


W  sin<p  sin0  COS0  -  U  COS0  COS0  -  U  cos<p  cos<p  sin0  sin0 
-  V  cos  <p  sin  9  sin  0  sin  0 

•  -Ucos(psin9sin0sin0- Vsin<j>sin<psin0sin0-V  cos0cos0  ■ 
-Wcos9cos0sin0 

-V  COS9  COS0  sin0  -  W  sin0  sin0  +  U  sin9  cos0  sin0 


A 

b  = 


1 

|U|sina 


V  COS0  +  W  COS9  sin  0 
■  Wsin9sin0-Ucos0 
-U  COS9  sin0  -  V  sin9  sinO 


xyt 


t  = 


sin  9  sin  0 
•  -cos9sin0 
COS0 


xyz 


tU 


cosa  =  -joj-  =  (Usin  9  sin0  -  V  cos9  sin  0 + W  cos0)/|U|. 


(14) 


THE  CURVATURE-BASED  COORDINATE  SYSTEM  (N,B,t) 

A  curved  line  in  space  can  be  used  to  define  three  directions:  tangent,  normal,  and 
binormal.  Consider  figure  6. 
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t  =  3 

N  =  pt'  p^O 

*  A  A 

B  =  TxN 

N 

Figure  6.  The  Curvature-Based  Coordinate  System 

A 

The  tangent  vector  T  in  this  coOTdinate  system  is  the  same  as  that  in  the  flow-based 
coordinate  system;  it  is  also  defined  as  3,  the  body-fixed  tangent  vector. 

The  noimal  vector  N  is  defined  by  the  curvature;  in  other  words,  by  the  way  T  changes. 
Locally,  the  curvamre  of  the  line  defines  a  plane,  the  plane  of  curvature,  f '  is  orthogonal  to  f , 
but  its  magnitude  is  equal  to  the  curvature,  not  to  unity.  Thus  the  reciprocal  of  the  curvature,  the 

A  A 

radius  of  curvature  p,  is  used  to  form  a  unit  vector  out  of  T'.  Finally,  the  binormal  vector  B  is 
fmmed  fiom  the  cross  product  of  f  and  N.  This  completes  our  orthonormal  set  ^N,  B,  f  j.  The 
N/B  plane  is  called  the  normal  plane,  the  B/T  plane  is  called  the  osculating  plane,  and  the  T/N 
plane  is  called  the  plane  of  curvature. 
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THE  DYNAMIC  EQUATIONS  ~  TRANSLATION 
Consider  a  differendal  secdon  of  cable  of  stretched  length  di,  as  shown  in  figure  7. 


Figure  7.  Force  Balance  on  Cable  Segment 


The  external  forces  acdng  on  the  cable  are  the  weight  per  unit  length  w  and  the  hydrodynamic  load 
per  unit  length  D .  The  internal  forces  acting  on  the  cable  segment  ate,  in  our  body-fixed 
coordinate  system  (l,  2,  §),  the  tension  f ,  and  two  shears,  %  and  V2.  Thus,  a  simple  force 
balance  yields 


Bt  31  3t 


di  -  wzdi  +  Dd/  =  ^df^. 


(15) 


In  equation  (15)  we  use  our  sign  convention  wherein  forces  are  considered  positive  when 
drawn  in  a  positive  direction  on  a  positive  face.  Both  w,  the  weight  per  unit  length,  and  |i,  the 
mass  per  unit  length,  arc  written  in  terms  of  the  stretched  length,  di.  We  can  write  everything  in 
terms  of  the  unstretched  length,  (Uq,  by  introducing  the  strain,  e,  and  the  stretch,  e: 
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e  = 

d/fl  d^o 


c  =  ^  =  e.l, 
d«o 


Thus,  we  can  introduce  the  unstretched  mass  and  weights  per  unit  length,  and  wq  , 


=  M 


Wo 

w  = 


and  rewrite  equation  (15)  as 


d^g  BIq  BIq 


-wo£  +  eD  =  iio® 


yf  +  V2  +  T'  -  Wo  2  +  eD  =  iio  u»  (21 

where  the  "  refers  to  differendadon  with  respect  to  /g,  and  the refers  to  differendadon  with 
respect  to  dme.  Equadon  (21)  is  decepdvely  simple  because  only  the  wg  term  and  the  u  term  are 
in  the  inerdal  coordinate  system  (x,  y,  z).  Written  explicidy  in  terms  of  the  body-fixed  and  flow- 
based  coordinate  systems,  equadon  (21)  becomes 

v/i  +  Vii' +  V2'2 +  rs  +  Ty 

(22: 

-Wg  z  +  eD„n  +  eD,t  =  lig(ux  +  vy  +  wz). 
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where  D„  is  the  ncrmal  component  of  the  drag  force  per  unit  length  and  is  the  tangential 
conqwnent,  both  to  be  specified  later. 

To  rewrite  equation  (22)  in  terms  of  only  the  fixed  reference  fiame  (x,  y,  z),  we  must  first 
replace  the  derivatives  of  the  body-fixed  unit  vectors  (l',  2\  by  their  values  given  in 

(AAA 

1, 2, 3,  n,  1 1  by  their  values 

given  in  equations  (4)  and  (14).  The  result  is  our  three  translational  dynamic  equations  of  motion: 

F,  =  m  a, 

T'(sin<psin0) 

+  T  (9'  COS9  sin0  +  0'  sin9  cos0) 

+  V/  (cos9  cosy  -  sin9  cos0  sintjr) 

+  Yi  (“9^  sin 9  cosy  -  v'  cos 9  siny  -  9'  COS9  cos0  siny 

•f  0' sin  9  sin  0  sin  y-y' sin  9  cos  0  cosy)  (23) 

+  V2  (-COS9  sin  y  -  sin  9  COS0  cosy) 

+  V2  (9'  sin 9  sin y  -  y'  cos 9  cosy  -  9'  c6s9  cos0  cosy 
-i-  0'  sin  9  sin  0  cosy  +  y'  sin  9  cos0  sin  y) 

-t*  —  f-U  cos^  9  sin^  0  -  U  cos^  0  -  V  cos9  sin9  sin^  0  W  sin9  cos0  sin0) 

fU|sina'  ' 

+  eD,  (sin  9  jin0) 

=  )ioU. 
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Fy  =  m  Cy 

T  (-cos(p  sin0) 

+  T  (<p'  sin<p  sin0  -  0'  cos<p  cos0) 

+  W{  (sin  <p  cos  V  +  cos<p  cos0  sin  y) 

+  Vi  (9'  cosq)  cosy  -  v'  sin<p  sin\jr  -  <p'  sinq)  cos0  sin  v 

-0'cos9sin0sinv  +  Y'cosq>cos0cosv)  (24) 

+  Yj  (-sin  9  sin  V  +  COS9  cos0  cosy) 

+  Vj  (-  9'  COS9  sin  V  -  v'  sin9  cosy  -  9'  sin9  cos0  cosy 
-  0' cos  9  sin  0  cos  V  -  v' cos  9  cos  0  sin  v) 

+  (-U  COS9  sin9  sin^ 0  -  V  sin^ 9  sin^ 0  -  V  cos^ 0  -  W  COS9  cos0  sin0) 

+  eDt  (-cos9sin0) 

=  m)  V. 


F,  =  m  a. 


rcos0 
+  T(-0'sin0) 

+  Vi(in0sinv) 

+  V|  (O'  COS0  sin  y  +  y'  sin  0  cosy)  (25) 

+  Vj  (sin  0  cosy) 

+  V2  (O'  cosO  cosy  -  y'  sin 0  sin  y) 

+  (u  sin9 COS0  sinO  -  V COS9  cosO  sinO -  W  sin^ o) 

+  eDt  (cos0) 

-Wo 

=  m)  w. 
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THE  DYNAMIC  EQUATIONS  --  ROTATION 


Consider  a  differential  section  of  cable  as  shown  in  figure  8. 


Figure  8.  Moment  Balance  on  Cable  Segment 

As  shown  in  figure  8,  the  cable  segment  is  subjected  to  the  internal  shearing  forces  and  tension 
(Vi.Vj.andf),  as  well  as  the  internal  bending  and  twisdng  moments  [Mj,  M2,  and  By 
taking  moments  around  the  center  of  the  segment,  we  need  never  introduce  the  weight  and  drag 
forces  (whose  torque  contribudons  are  of  a  higher  OTder).  In  the  body-fixed^rincipal-axis 
coemlinate  system,  the  dynamic  equations  of  rotation  are 


ti  =  Ii  ©1  +  (I3  - 12)  ©2 

X2  =  I2  ©2  +  (^1  "•  I3)  ®3  (26) 

T3  =  I3  ©3  +  (I2  “  Ii)  ©I  ©2 , 

where  is  the  torque  about  the  ith  axis;  Ij,  the  principal  moment  of  inertia  about  the  ith  axis;  and 
©i,  the  angular  speed  about  the  ith  axis. 
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From  figure  8,  the  net  tOTque  on  our  differential  segment  is 


dx  =  — J-df +  -:^df +  ^r^df  + V,df  2- V2df  1. 
dl  didi 


(27) 


By  introducing  the  unstretched  length  [equation  (17)],  equation  (27)  is  changed  to 


dx  =  M[d/o  +  Mid/o  +  M[dfo  +  Vjed/o2-V2edfol 


(28) 


V  =  M[  +  Mi  +  M(  +  Vie2-V2ei  (29) 

or 

x'  =  M[i  +  Mir  +  Mi2  +  M2r+M;3  +  M,3'  +  Vie2-V2ei.  (30) 

The  three  principal  moments  of  inertia  of  our  segment  of  length  dt  are; 
dli=dl2=^lid/r2 

dl3=^Hd/r2.  (31) 

Again,  we  introduce  the  unstretched  lengths 
IJ  =  =  Vo, 

l5=2IS,and  (32) 

lo  ■  iUor^- 
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In  the  body-fixed/principal-axis  coordinate  system,  the  angular  velocity  is  written  in  terms 
of  the  Euler’s  angles  as  (see  reference  1): 


(0  =  ( 


<p  sin0  sin\|r  +  9  cosy 
(p  sin6 cosy  -  Osin y 
(pcos0  +  y 


123 


(33) 


Hnally,  using  equation  (33)  for  the  angular  velocities,  equation  (32)  for  the  moments  of 
inertia,  equation  (30)  for  the  ttxques,  and  equation  (13)  for  die  derivatives  of  the  unit  vecton,  the 
governing  dynamic  equations  of  rotation  (26)  can  be  rewritten  as; 

I  M(  +  M2  (-<p'  COS0  -  y')  +  Mj  (9'  sinO  cosy  -  0'  siny)  -  e  V2 

-Ifl  |9sin0siny  + 0  cosy +  <p^  COS0 sin 0 cosy  (34) 

+  2  <p  y  sin0  cosy  -  2  0  y  siny]  =  0 

i 

I  Ml  (9'  COS0  +  y')  +  M2  +  M{  (-  9'  sin 0  sin  y  -  0'  cosy)  +  e  Vi 
I  -Ifl  |9sinOcosy-0siny-9^  cosOsinOsiny  (35) 

I 

-  29y  sin0  siny  -  20y  cosy|  =  0 


Ml  (-9'sin0cosy +  0'siny)  +  M2  (9'sin0siny +  0'cosy) 

+  MJ-2Io{9Cos0-90sin0  +  yJ  =  O.  (36) 
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CONSTITUTIVE  AND  KINEMATIC  RELATIONS 

In  the  second  and  third  secdons,  we  developed  six  dynamic  equations  in  13  unknowns: 
three  moments  (Mj,  M2.  Mj,  three  forces  (Vj,  V2,  T),  three  velocities  (u,  v,  w),  three  angles 
(9. 6,  y),  and  the  stretch  e.  In  this  secdon  we  will  develop  the  seven  equadons  needed  to  close 
the  system. 


CONSTITUTIVE  EQUATION  FOR  THE  TWISTING  MOMENT 

According  to  the  linear  theory  of  elasdcity,  the  twisting  moment  and  die  twist  per  unit 
length  are  linearly  related  through  the  shear  modulus  and  the  polar  cross-secdonal  moment: 

M,  -  (37) 

If  we  assume  that  the  caUe  is  not  torque-balanced,  then  we  can  write  that  the  twisdng  moment  is 
also  linearly  related  to  the  strain  e  ~  e  - 1: 


M,  =  A4  y'  +  A3(e-I). 


(38) 


CONSTITUTIVE  EQUATION  FOR  THE  TENSION 

Again,  according  to  the  linear  theory  of  elasdcity,  the  tension  and  the  strain  are  linearly 
related  through  Young's  Modulus  and  die  cross-secdonal  area: 

T  =  EAe. 
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Again,  if  we  assume  that  the  cable  is  not  torque-balanced,  then  we  can  write  that  the  tension  is  also 
linearly  related  to  the  twist: 

T  =  A,  (e-l)-»-A2  V'.  (39) 

CONSTITUTIVE  EQUATIONS  FOR  THE  BENDING  MOMENTS 

According  to  pure  bending  thecay,  which  we  will  use  despite  the  fact  that  shears  are 
present,  the  bending  moment  is  linearly  related  to  the  bending  rigidity  (E  I),  and  inversely  related 
to  the  radius  of  curvature  (p).  Figure  9  shows  this  reladonship  in  terms  of  the  ^N,  B,  f  j 
coordinate  system,  the  coordinate  system  based  on  the  local  curvature,  developed  in  the  first 
section. 


Mb 


Figure  9,  Pure  Bending 
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T7 


Because  the  bending  moment  is  defined  as  being  in  the  binormal  direction,  we  can  write, 
from  figures  6  and  9, 

Mb  = -^B  =  ElBTx^t  (40) 

or 


cMb  =  ElBtxr. 


Using  equation  (13)  for  f '  =  33/3fo ,  we  can  write 


txT'  as 


■  (-  <p'  sind  siny  -  9'  cosy)] 
9'  sinO  cosy  -  6'  siny 
0 


123 


(41) 


(42) 


or 

cMb  =  EIb 

A 

where  die  cross-sectional  bending  moment  Ib  is  fcH*  bending  around  the  B  axis  and  is,  for  a 
homogeneous  cable,  icR^/4.  Thus,  the  constitutive  equadons  for  the  two  bending  moments  are 


9'sin0siny +  9  cosy 
9'  sin9cosy-9'  siny 
0 


(43) 


eMj  =  EIb  (9' sin  9  siny +  0' cosy) 


(44) 


and 
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eM2  =  Elg  (9'sin0cosv-0'siny). 


(45) 


KINEMATIC  RELATIONS 


The  (x,  y,  z)  coordinates  of  a  point  on  the  cable,  and  the  (9, 0,  v)  angles  at  that  point,  are 
independent  However,  the  way  in  which  the  coordinates  change  along  the  cable  is  determined  by 
die  orientation  of  the  cable,  by  the  angles.  Specifically,  the  cable  coordinates  change  in  the 
direction  of  the  tangent  vector 


9T 

di 


(46) 


which,  written  out  by  components,  gives  us  three  kinemadc  equations: 
x'  =  esin9sin0, 
y'  =  -ecos9sin0. 


(47) 

(48) 
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SUMMARY 


In  this  report  we  developed  a  complete  set  of  field  equations  describing  the  dynamic 
behavior  of  a  towed  cable  and/m*  array.  The  cable  has  bending  stiffness,  tOTsional  sdffness, 
extensional  stiffness,  inertia,  and  viscous  drag.  The  viscous  drag  (as  well  as  added  mass)  is  left  as 
a  generic  term  to  be  defined  by  the  users  according  to  their  specific  needs. 


The  ctmplete  set  of  field  equadons  cmisists  of  13  nonlinear  partial  differential  equations 
and  is  repeated  here  for  conciseness; 


L 


F* 


m  a. 


T  (sin  (p  sin  6) 

+  T  (<p'  cos  jp  sin  0  +  0'  sin  <p  cos0) 

+  V{  (cos<p  cos -  sin 9  cos 0  sin  v) 

+  Vj  (-<p'  sin9  cosy  -  v'  cos<p  sin y  -  <p'  coscp  cos0  sin y 
+  9'  sintp  sin0  siny  -  y'  sin9  cos0  cosy) 

+  Vj  (-COS9  sin y  -  sin  9  cos0  cosy) 

+  V2  (9'  sin9  siny  -  y'  cosy  cosy  -  9'  cosy  cos0  cosy 
+  0'  siny  sin0  cosy  +  y'  siny  cos0  siny) 

+  ■  f-U  cos^  9  sin^0  -  U  cos^0  -  V  cosy  siny  sin^ 

(Ulsino' 

+  eD,  (sin9sin0) 


(50) 


0  +  W  siny  cos0  sin0) 
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2. 


Fy  =  m  ay 


T' (-cos<p  sin0) 

+  T  (9'  sin  9  sin  0  -  0'  C0S9  cos0) 

+  y{  (sin9  cosy  +  COS9  cos0  sin  9) 

+  Vj  (9'  cos  9  cos  9  -  9'  sin  9  sin  9  -  9'  sin  9  cos  0  sin  9 

-  0' cos  9  sin  0  sin  9  +  9' cos  9  cos  0  cos  9) 

+  yj  (-sin  9  sin  9  +  COS9  COS0  cos  9) 

+  V2  (-  9'  COS9  sin  9  —  9'  sin  9  COS9  —  9'  sin  9  cos0  C0S9 

-  0' cos  9  sin  0  cos  9  -  9' cos  9  cos  0  sin  9) 

+  (-U  COS9  sin9  sin^0  -  V  sin^9  sin^0  -  V  cos^0  -  W  cos^) 

jU|sina' 

+  cD,  (-cos9sin0) 


(51) 


cos0sin0) 


=  (lio+li.)v- 


3.  E,  =  m 

rcos0 
+  T(-0'sin0) 

+  y((sin0sin9) 

+  Vi  (0' COS0  sin9  +  9' sin0  COS9)  (52) 

+  V2  (sin  0  cos  9) 

+  Yz  (0'  COS0  C0S9  -  9'  sin 0  sin 9) 

+  (U  sin9 COS0 sin0  -  V COS9 cos0 sin0 -  W sin^0) 

+  eD,  (cos0) 

-Wo 

=  (1^0+  lit)w- 

4.  MJ  +  M2  (-9' COS0 -  9')  +  M,  (9' sin0 C0S9  - 0' sin9) - c V2 

-Iq  {9sin0sin9  +  0cos9  +  9^  cos0sin0cos9  (53) 

+  299  sin0  C0S9  -  2  0  9  sin9|  =  0 . 
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8. 


10. 


13. 


M|  (<p'  cos9  +  v')  +  +  Mj  (-  (p'  sin  0  sin  v  -  0'  cosy)  +  e  Vj 

-  lo  sin  0  cos  V  -  0  sin  V  -  <p"  COS0  sin  0  sin  v 

-  2<p  V  sin0  sin  V  -  20\ir  cosvj  =  0 . 

Ml  (-  <p' sin0  cosy  +  0' sin y)  +  M2  (<p' sin0  sin y  +  0' cosy) 
+  M(-2Io  |9cos0-«j)0sin0  + \jir|  =  0. 

M,  =  A4  v'  +  A3(c-1). 

T  =  AiCe-O  +  AaV'. 


9.  eMi  =  Elg  (<p'sin0sinv  +  0'cos\/). 


eM2  =  EIb  (9'sin0cosv-0'sinv). 


IL  u'  =  ~(esin9sin0). 


12.  v'  =  -  — (ecos9sin0). 


w'  =  ^(ecos0). 


(54) 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


The  13  unknowns  are  the  3  forces  (Vi,  V2,  T),  3  moments  (M|,  M2,  Mj,  3  velocities 
(u,  V,  w),  3  angles  (9, 9,  y),  and  the  stretch  e.  Written  as  such,  the  13  equations  (50-62)  are  first 
Older  in  space  and  second  order  in  time.  By  introducing  the  3  angular  velocities. 


X  =  <P. 


(63) 


24 


/ 


/ 


TR  10,231 


1  =  6,  and 


©  =  V, 


we  would  increase  our  system  to  16  equadons  in  16  unknowns  and  become  first  order  in  space  and 


Hnally,  because  these  equadons  will  potentially  be  used  to  study  the  higher  order  effects  of 
rigidity  and  torsion,  we  can  simplify  these  equadons  for  comparison  and  reduce  them  to  the  more 
standard  six  unknown  systems  written  here  in  equadons  (66)  through  (71): 

(sin  <p  sin  6) 

+  T  (<p'  cos<p  sin  0  +  0'  sin  <p  cos0)  (66) 

+  ^  (*"U  cos^  <p  sin^  0  -  U  cos*  0  -  V  cos <p  sin 9sin*  6  +  W  sin <p  cos0  sin 0) 


+  D,  sinq)  sin0  =  (|io  +  u 


T' (-cos  9  sin  0) 

+  T  (9'  sin  9  sin  0  -  0'  cos9  cos0)  (( 

+  ^  (-U COS9 sin9 sin*0 -  V sin*9 sin*0 -  V cos*0 -  W COS9 cos0 sin0j 


+  D,  (-cos9sin0)  =  (^o+  ji,)v 


rcos0  +  T(-0'sin0) 

+ (u  sm9  cos6  sin0  -  V  cos9  cos0  sin0  -  W  sin*0j 


+  D,  COS0-WO  =(lio  + 


u  =  ^  (sin  9  sin  0) 
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- (cos  (D  sin  6) 

dt 


(70) 


—  (cos0) . 
dt 


(71) 


The  soludons  to  these  equations  will  be  the  subject  of  a  subsequent  report 


I 

I 

I 

I 

I 

I 

I 

I 

I 
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